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DOES DISK LOCKING SOLVE THE STELLAR ANGULAR MOMENTUM PROBLEM? 

Sean Matt^-^ and Ralph E. Pudritz"' 
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ABSTRACT 

We critically examine the theory of disk locking, which assumes that the angular momentum de- 
posited on an accreting protostar is exactly removed by torques carried along magnetic field lines 
connecting the star to the disk. In this letter, we consider that the differential rotation between the 
star and disk naturally leads to an opening (i.e., disconnecting) of the magnetic field between the two. 
We find that this significantly reduces the spin-down torque on the star by the disk. Thus, disk-locking 
cannot account for the slow rotation (^^10% of breakup speed) observed in several systems and for 
which the model was originally developed. 

Subject headings: accretion, accretion disks — MHD — stars: formation — stars: magnetic fields — 
stars: pre-main-sequence — stars: rotation 
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1. INTRODUCTION 

We now know that most protostars undergo a phase in 
which they accrete mass from a disk and that a stellar 
magnetic field dominates this process near the star. How- 
ever, the transport of angular momentum, while fairly 
well-understood within the disk, remains a significant 
challenge for mode ls of stellar spin (for a review, see 
iBodenheimeU 119951) . Classical T Tauri stars (CTTS's) 
comprise a class of accreting protostars for which there is 
an abundance of observational data. For typical param- 
eters in these systems, the angular momentum deposited 
by accretion can spin the star up to breakup speed in 
~ 10^ years, assuming the star hadn't already formed 
at breakup speed. Since the accretion lifetime is often 
greater than 10^ yr, the stars must somehow rid them- 
selves of this excess angular momentum. Furthermore, a 
significant number of CTTS's (the so called "slow rota- 
tors") have spin rates as slow as ^10% of breakup speed 
(e.g., [Bouvier et al. 1993J, requiring significant torques 
t o maintain. 

iKonigM l|19911 hereafter K91) proposed the currently 
accepted view that the slow rotators could be ex- 
plained with a "disk-locking" (DL) mechanism (adapted 
from Ghosh & Lamb 1979), for which the star needs 
a kilogauss-strength dipole magnetic field. Indeed, 
CTT S's are now known to possess kC-stre ngth fields 
Ce.g.. lJohns-Krull. Valenti. fc Koreskolll999bj) . and mag- 
netic star-disk interaction models have also explained 
several observed spectral featur es of CTTS 's (e.g., 
K91 ). However, observations by iStassun et al.l (|199§) 
and I.Tohns-Kr^iU et alJ |1999a^ and theoretical consid- 
erations of iSafierl l)1998t) have called the standard DL 
scenario into question. According to DL theory, mag- 
netic field lines causally connect the star to the disk 
(acting as "lever arms") and carry torques that oppose 
and balance the angular momentum deposited by ac- 
cret ion. Although not usually consid ered in DL the- 
ory, lUzdenskv. Konigl. fc L itwinI l)2002l hereafter UKL) 
demonstrated that differential rotation between the star 
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Fig. 1. — Magnetic star-disk interaction. 



and disk leads to an opening of the field, drastically re- 
ducing the magnetic connection between the two. In this 
letter, we consider the consequences of this reduced con- 
nection on DL. We find that the magnetic torque is signif- 
icantly reduced, and show that the DL model therefore 
cannot account for the angular momentum loss of the 
slow rotators. 

2. TORQUES IN THE STAR-DISK INTERACTION 

In order to discuss the disk-locked state, in which mag- 
netic torques on the star exactly balance the angular mo- 
mentum deposited by accreting material, we must first 
discuss the general case for an arbitrary state of the 
star-disk system but with the proper inclusion of field 
line openin g. To this e nd, we mainly follow the for- 
mulation of lArmitage fc Clarke (1996, hereafter AC96), 
and a more detailed pre s entati on of our model will be 
given in lMatt fc Pudritd l|2004D . According to the typi- 
cal assumptions, the central star contains a rotation-axis- 
aligned dipolar magnetic field. This field is anchored in 
the stellar surface and also threads a thin, Keplerian ac- 
cretion disk, in which the kinetic energy of orbiting gas 
is much greater than the magnetic field energy. The re- 
gion above the star and disk contains low-density gas 
and is dominated by the magnetic field. The system is 
assumed to exist in a steady-state, where the accretion 
rate Ma is constant in time and at all radii in the disk. 
Figurenillustrates the basic idea and identifies the loca- 
tion where the disk is truncated {Rt) and the corotation 
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radius (i?co), where the Keplerian angular rotation rate 
equals that of the star. The region in the disk threaded 
by closed stellar field lines is delimited by i?out (the usual 
assumption is i?out ^ -Rco)- 

The assumption of a steady-state Ma requires that the 
net angular momentum carried away from each annulus 
of the disk (of width 5r and vertically integrated) equals 

6t, = 0.5 AU^GmI (i) 

where A/, is the stellar mass. If the torque anywhere 
differs from this, the disk will restructure itself on an 
dynamical timescale, reestablishing the steady-state. 

The star and disk rotate at different angular speeds, 
except at the singular location i?co- Thus the mag- 
netic field becomes twisted azimuthally by differential 
rotation, and magnetic forces act to restore the dipole 
configuration — conveying torques between the star and 
disk. The magnetic torque on the star exerted by a bun- 
dle of field lines threading an annulus in the disk, can be 
written as (AC96) 

5T^,^^^?r~Hr, -i = B^/B,. (2) 

This differential torque depends only on the strength of 
the dipole moment /i and on 7, which is the "twist" (or 
pitch angle) of the field. The sign of the torque, here 
and in all equations, is relative to the star, so a positive 
torque spins the star up and thus spins the disk down. 

The magnetic field cannot be perfectly frozen-in to 
the gas of the disk. Instead, it diffuses or reconnects 
through the disk azimuthally at some speed, Wd ~ 7 ?7t//i 
(UKL), where rjt is the magnetic diffusivity (we use the 
subscript "t" to suggest a turbulent value) and h is the 
local, vertical scale height of the disk. For a standard 
a-disk (^hakura & Sunyaev 1973), Wd — /3t^kcp7, where 
Ukcp is the orbital speed, (3 = {a/Pt){h/r), and Pt is the 
turbulent Prandtl number (= turbulent viscosity divided 
by ?7t). Since both h/r and a/Pt are likely to have weak 
(though unknown) dependences on r, we assume /3 is con- 
stant. The radial dependence of the magnetic torque is 
mainly dominated by the quick falloff of the dipole mag- 
netic field {r~^), and so a small radial dependence of (3 
will not much affect our results (AC96). 

In general, /3 is a simple scale factor that compares 
to Wkop- It measures the coupling of the magnetic field 
to the disk: /3 <C 1 corresponds to strong coupling and 
/3 » 1 to weak. The value of f3 is unknown (AC96 used 
/3 = 1). Typical a-disk parameters give a limit of /3 < 1 
and a likely value of a few orders of magnitude lower. 
For reasonable fiducial parameters, /3 — 10~^(77t/10^^ 
cm^ s~^)(i?0//i)(lOO km s~^/ukcp)- However, given the 
uncertainties, we keep /3 as a free parameter. 

Where the field connects the star and disk, the mag- 
netic twist (and so Vd) will increase, until the field can 
slip through the disk at a rate equal to the differential 
rotation rate. Thus, we expect a steady-state in which 



7 = r 



(3) 



Note that 7(i?co) — 0, since the differential rotation is 
zero there. For r > i?co, the twist increases to infinity. 

Several recent theoretical and numerical studies of the 
star-disk interaction (s ee UK L and references therein; 
Lvnden-Bell &: Boilvl 119941: lAgapitou fc Papaloizoul 
20001) have shown with certainty that, as a dipole field 
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Fig. 2. — Differential torques in the disk midplane predicted for 
/3 = 1 and 7c — > 00, for BP Tau (see text). The system is shown 
in the equilibrium state, where the net torque on the star is zero, 
requiring a stellar rotation period of 7.5 d (so iico ~ 6.4i?*). 



is twisted this way, the magnetic torque reaches a 
maximum for finite 7, but then reduces back to zero for 
larger values. This occurs because, when the magnetic 
field is twisted enough so that the magnetic pressure 
associated with overcomes the poloidal field lines, the 
latter will "inflate," opening to infinity at mid latitudes. 
The star and disk then become causally disconnected 
because open field lines cannot convey torques between 
the two. The critical twist at which this happens is very 
nearly 7 = 7^ « 1 (e.g., UKL). 

The opening of field lines by this process can be added 
to the theory. For simplicity, we assume that the star is 
connected to the disk, as described above, at all locations 
where the twist is less than the critical value 7c. Where 
7 7c, the disk and star are disconnected, and so the 
differential magnetic torque there is zero. In other words, 
this process determines the outer radius of closed stellar 
field lines, i?out- From equation 13 we find that i?out = 
(l-|-/97c)^^"^i?co- AC96 assumed that the star is connected 
to the disk over a very large radial extent (i?out 00), 
which is equivalent to the assumption that 7c — > cxd. 

A combination of equations |5] and 13 reveals the ra- 
dial dependence of the differential magnetic torque (St^). 
Figure [21 shows all differential torques (per Sr) in the 
disk, as a function of radius (normalized to i?co), assum- 
ing (3 = 1 and 7c — > cxo (i.e., the AC96 solution). For 
the Figure, we adopt representative system parameters 
from the well-studied CTTS BP Tau: = 3 x lO^fM^ 
yr"^ = 2i?Q, and = O.SMrr^ IGullbring et alJ 
11991 . and = 2 kG (|.Tohns-Krul] TTalTll 999H) . The 
stellar rotation period is 7.5 days (see ®. 



The lines in Figure |21 labelled "accretion" and "mag- 
netic" represent the differential accretion and magnetic 
torques, (Sra and (STm, from equations ^ a-nd [3 respec- 
tively. The steady-state condition requires the disk to 
structure itself such that the net differential torque ev- 
erywhere equals St^,. So when external magnetic torques 
(Jthi) act on the disk, there are torques internal to 
the disk (e.g., via th e magnetorotational instability, 
iBalbus & H awlcv 1991) to counteract them and always 
provide a differential torque defined by Sti = (5Ta — (5Tm, 
represented by the dashed line in Figure From the 
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figure, it is evident that is strongest near the star, 
where the magnetic field is strong, and Jrm acts to spin 
up the star inside Rco- At i?co, goes to zero, since 
the field is not twisted there (eq. Beyond Rco, the 
twist increases, and so St^ becomes stronger, now act- 
ing to spin down the star. Since the dipole field strength 
decreases faster than the increase of 7, Jrm reaches a 
maximum (in absolute value) and then approaches zero 
as r ^ 00. 

At the location where St^ = St^ (and St^ = 0), the 
torque from the stellar magnetic field is all that is necce- 
sary to provide Mg_. Thus, all of the specific angular 
momentum of the disk material at that location will end 
up on the star. This defines the radius i?t (dotted line in 
Fig. 121), where the disk is truncated, and from where ac- 
cretion will be magnetically channelled along field lines 
onto the star (K91; not shown in Fig. So the net 
torque on the star from the accretion of disk material, 
Ta, is obtained by integrating equation ^ from i?t to the 
surface of the star. Similarly, the net magnetic torque. 
Tin, is obtained by integrating equation |5] over the entire 
magnetically connected region of the disk, from Rt to 
i?out, which gives 

T,„ = iJi\iP)~^R-^' [2 (1 + /37c)"' - (1 + /37c)"' 

~2{RcolRtf'^ + [Rco/Rtf] ■ (4) 

This is exactly the solution found by AC96 for the special 
case that P — 1 and 7c — > cxd. Our formulation allows 
for an arbitrary value of the diffusion parameter /3 and 
includes the effects of field line opening via twisting to a 
critical value of 7c . 

We determine the dependence of the spin-down torque 
on P by setting Rt = Rco in equation ^ and adopting 
7c = 1 (as justified by, e.g., UKL). We find that Tm cx /3 
in the strong coupling limit (/3 <C 1), has a maximal value 
for l3 — 1, and oc /3~^ in the case of weak coupling 
(/? ^ 1). This behavior can be understood as a compe- 
tition between two different effects: one is that, when /3 
is small, the size of the magnetically connected region in 
the disk (i?out) is small, reducing t,„; secondly, when /3 
is large, the field is less twisted (eq.EJ at each radius, re- 
ducing St^ . These two effects conspire to give a maximal 
value of Tm, for the critical value of P — 1, which thus 
represents the "best case" for DL theory. Even for this 
"best case," the proper treatment of field line opening 
(7c = 1) gives a total spin-down torque that is a factor 
of four times less than if 7c —^ 00. 

To illustrate the effect of field line opening. Figure 13 
shows the differential torques in the disk, for the same 
parameters as Figure |21 except that now 7c = 1, and the 
spin period is 4.1 d (see The magnetic field is now 
open for r > i?out ~ l-6i?coj and Srm/dr is zero there 
(by assumption). It is clear that the integrated torque 
Tin will be less than if the magnetic field were everywhere 
closed (i.e., for 7c — > 00). 

3. THE DISK-LOCKED STATE 

For any given values of M,, i?,, _B*, Afa, and the stel- 
lar rotation period, the star-disk interaction theory pre- 
sented in section 121 gives the net torque (ra + t„^) on the 
star. The system is stable in that, for fast rotation, the 
net torque spins the star down, and for slow rotation, 
the star spins up. For typical CTTS parameters, the net 
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Fig. 3. — Same as Figure I2I except that 7c = 1, so the star is 
magnetically disconnected from the disk beyond Rout ~ 1.6i?co, 
and the equilibrium rotation period is 4.1 d (-Rco ~ 4.3/?*). 



torque will spin the star up or down in ^ 10^ yr (e.g., 
AC96). So one expects that most systems older than 
this will exist in an equilibrium state, in which the net 
torque on the star is zero, and in which the system is 
"disk locked" with a period Tcq. 

The equilibrium spin state is determined by combin- 
ing the condition of net zero torque (ra -I- Tm = 0) with 
the definition of Rt (i.e., where (5ra = (^Tm). One finds 
that, in the DL state, there is a single predicted value of 
the truncation radius, i?^'', that depends on the quantity 
/37c, but not on any other system parameters. If /?7c ap- 
proaches zero, Rt"^ approaches Rco- The ratio Rt'^/Rco 
monotonically decreases for increasing Pjc, but it never 
becomes smaller than 0.91, as /37c 00. 

We calculate that the equilibrium angular spin rate of 
the star is given by 



= C(/3, 7c) Mf/' (GM.f B^'^'^R-^^'^ 



(5) 



where C(/3, 7c) is a dimensionless function that depends 
only /3 and 7c. 

Figures (21 and O show the torques for stars that are 
in their equilibrium spin states, Tcq = 7.5 and 4.1 d, re- 
spectively. A comparison of these two Figures reveals the 
effect of field line opening (see When this occurs, the 
net spin-down torque is decreased, so the star must spin 
faster for it to be in equilibrium. A faster spin reduces 
Rco, so the torques come from closer to the star, where 
the dipole field is stronger. Thus, a faster spin makes up 
for the decreased size of the magnetically connected re- 
gion. Also, the truncation radius is nearer the corotation 
radius in Figure than in|21(i?t'' ~ 0.97i?co, compared 
to « 0.91i?co). 

Equation has the same dependences on system pa- 
rameters as_sevBraJ_^L models in the literature (e.g., 
K91 and'Ostriker fc ShdlT995^ . but the factor C varies 
slightly from model to mo del. For example, both K91 
and Ostriker fc Shil l|T99l use C 1.1, while the AC96 
value (i.e., for /3 = 1 and 7c — > 00) is C « 1.6. Our 
formulation of the problem includes the effects of field 
line opening and determines the dependence of the func- 
tion C on the diffusion parameter /?. The solid line of 
Figure ^ reveals this dependence, for 7c = 1. For com- 
parison, if the field remains closed for arbitrary twist 
values (7c ^ 00), C(/3) « 1.6/3^/7 (dotted fine). The 
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Fig. 4. — Spin-rate factor C of equation |21 as a function of 
log(/3). The solid line represents the model with field line opening 
(7c = 1). For comparison, the dotted line represents 7c — » oo, 
while the dashed line is from K91. 



K91 value of C is indicated by a dashed line. The effect 
of field line opening is significant. For the "best case" 
(/3 = 1), the predicted spin rate is 1.8 times faster than 
for 7c — > oo and 2.6 times faster than predicted by K91. 
For a more reasonable value of /? = 10~^, the predicted 
spin is more than an order of magnitude faster than all 
other DL models, in which the field is typically assumed 
to remain closed while largely twisted. 

To date, DL models have had some success in explain- 
ing the spin rates of slow rotators. For the example 
case of BP Tau (see 33), when one assumes /3 = 1 and 
7c 00, equation O predicts Tcq — 7.5 d, corresponding 
to 6% of breakup speed and agreeing with the predic- 
tions of most models in the literature. This predicted 
pe riod is remarkably similar to the observed value of 7.6 
d ijVrba et alJll986|) . However, when one properly in- 
cludes the effect of field line opening (for 7c — 1), the 
"best case" (/3 = 1) model predicts Tcq — 4.1 d, and 
the spin-up time is less than 10^ yr — sign ificantly shorte r 
than BP Tau's age of 6 x 10^ yr ( GuUbriJZiiaOIIMi) • 
Further, for f3 = 0.01, Tcq = 1.0 d, with a spin-up time 
of 5 X 10^ yr! 

It is also important to note that we h ave thus far used 
parameters for BP Tau determined by iGullbring et al.l 
}i998) and l.Tohns-Krull et all l|1999bt). If instead we use 
the higher derived bv'Harti gan et al"! l)1995(l of 1.6 x 

10^'^ M0 yr^^, even the "standard" (/3 = 1, 7c — > oo) 
model predicts Tcq = 3.6 d, with a spin-up time of 2 x 



10* yr. Alternatively, if we use the 3a upper limit on 
the dip ole field strength of 200 G from lJohns-Krull et alJ 
lll999aD . even the "standard" theory gives Tcq = 1.0 d. 

We must conclude that either BP Tau and other slow 
rotators are rapidly spinning up, or the DL picture, as 
an explanation of stellar spin, is incomplete. The former 
requires (e.g.) a recent, dramatic increase in Ma — this 
seems unlikely, especially for BP Tau, for which the ac- 
cretion disk is probably in the process of clearing out 
ijDutrev et al. 2003). In order for BP Tau to currently 
be in an equilibrium spin state, there must be significant 
spin-down torques on the star other than those carried 
along field lines connecting the star to the disk. 

4. SUMMARY OF PROBLEMS WITH DL 

We have shown that a large portion of the magnetic 
field connecting the star to the disk will open up due to 
the differential rotation between the two, resulting in a 
spin-down torque on the star by the disk that is much 
less than if the field is assumed to remain closed. The 
predicted disk-locked spin rate is therefore much faster, 
so that the DL scenario cannot explain the angular mo- 
mentum loss of the slow rotators. 

In addition, there are at least three, completely inde- 
pende nt issues raised by other authors: 1) Stassun et alJ 
(1999) found no correlation between accretion parame- 
ters and spin rates of TTS in Orion. 2) CTTS's appar- 
ently do not have strong dipole fields (e.g.,[^ficr 199^ 
iJohns-KruU et al.1ll999a|) . which are required for DL at 
slow spin rates. 3) Stellar winds are expected t o open 
field l ines that would otherwise connect to the disk ijSafieil 
Il998j) . A disk wind could have a similar effect. 

We conclude that, in order for accreting protostars to 
spin as slowly as 10% of breakup speed, there must be 
spin-down torques acting on the star other than those 
carried by magnetic field lines connecting the star to the 
disk. The presence of open stellar field lines implies that 
excess ang ular momentum may be carried by a stellar 
wind (e.g.. iTout fc PringlellT99l . but this remains an 
open question. We are investigating the role of stellar 
winds, since then Tcq will depend mostly upon the stellar 
wind parameters, thereby explaining the lack of correla- 
tion between spin periods and accretion parameters. 
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